Transfer-matrix description of heterostructures involving superconductors and 

ferromagnets 
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Based on the technique of quasiclassical Green's functions, we construct a theoretical framework 
for describing heterostructures consisting of superconductors and/or spin-polarized materials. The 
necessary boundary conditions at the interfaces separating different metals are formulated in terms 
of hopping amplitudes in a t-matrix approximation. The theory is applicable for an interface with 
arbitrary transmission and exhibiting scattering with arbitrary spin dependence. Also, it can be 
used in describing both ballistic and diffusive systems. We establish the connection between the 
standard scattering-matrix approach and the existing boundary conditions, and demonstrate the 
advantages offered by the t-matrix description. 

PACS numbers: 74.20.-z,73.40.-c,74.50.-|-r 



I. INTRODUCTION 



Low-energy (below the superconducting energy gap) 
electron transport through contacts between a supercon- 
ductor and a normal metal can be understood in terms 
of Andreev reflectioni In this process, an incident elec- 
tron from the normal side can enter the superconductor 
by pairing with another electron with the opposite spin, 
leaving a reflected hole in the normal metal. The phase- 
coherent nature of this process results in superconducting 
correlations being induced in the normal-metal side, re- 
ferred to as the proximity effect. The important feature 
of Andreev reflection is that, with singlet superconduc- 
tors, it involves both spin bands in the normal metal. 
Therefore, the above simple picture has to be modified 
when the normal metal is replaced by a ferromagnet with 
two different Fermi surfaces for the two spins, result- 
ing in new and interesting physical phenomena. In re- 
cent years, interplay between superconductivity and fer- 
romagnetism has attracted considerable theoreticalSiSi^ 
and experimentalSi^iS attention - both out of fundamen- 
tal scientific interest and in view of the possibility of novel 
applications and devices. One important consequence of 
the spin splitting between the two bands in the ferromag- 
net is that the phase coherence between the particle-hole 
pair in the clean (dirty) limit is destroyed over a charac- 
teristic distance of Vf /h (y^D/h), where Vf is the Fermi 
velocity, D the diffusion constant, and h an effective ex- 
change energy which describes the spin splitting. Since 
this distance is typically very short (of the order of atomic 
distances), the superconducting correlations induced to 
the ferromagnetic material are expected to be confined to 
the immediate vicinity of the separating interface. This 
raises the question whether, for strong ferromagnets, a 
mechanism of a different type takes over and dominates 
the physics of superconductor/ferromagnet (S/F) con- 
tacts. One such mechanism, recently under active inves- 
tigation, is the inducement of spin-triplet correlations: 
namely, the exchange field only affects correlations of 



singlet type, i.e between particles and holes in opposite 
spin bands. In fact, equal-spin triplet correlations are 
expected to be created by proximity to a ferromagnet 
due to the breaking of spin-rotational symmetry. The 
desire to formulate a theory capable of understanding 
the detailed nature and the conditions for the formation 
of these correlations and the corresponding anomalous 
proximity effect has given the initial motivation for this 
work. 

Problems related to superconducting proximity ef- 
fect with spin-dependent interfacial scattering are of 
spatially inhomogeneous nature and, as such, they 
can only be studied with specialized theoretical tools. 
Such a tool is provided by the quasiclassical theory 
of superconductivityiiSiii This theory is applicable for 
weakly perturbed superconductors (characteristic length 
scale of perturbations much larger than Fermi wave 
length and characteristic frequencies much less than 
Fermi energy) and can be used in both equilibrium 
and nonequilibrium situations. It describes quasiparti- 
cles with momenta on the Fermi surface moving along 
straight classical trajectories, the direction of which is 
given by the corresponding Fermi velocity. A ferromag- 
netic metal has different Fermi surfaces and, correspond- 
ingly, different sets of trajectories for the two spin ori- 
entations. In this case, the quasiclassical theory can be 
used to describe two limiting cases: i) weak ferromag- 
netism, where the energy splitting of the two Fermi sur- 
faces and the associated deviation of the Fermi velocities 
is so small that the two spin trajectories with the same 
momentum direction are fully coherent, and ii) strong 
ferromagnetism, where the splitting is so large that the 
coherence is lost completely. While the former limit- 
ing case has been exhaustively studied in the literature, 
the latter has only recently received attention,-'^^- Here 
we present a theoretical study of the latter possibility. 
Even in the absence of conventional Andreev reflection 
processes (which would require coherence between par- 
ticles and holes in opposite spin bands), interesting and 
nontrivial physics emerges due to spin-active interfacial 
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scattering. Additional motivation has been provided by 
the growing interest in a new class of materials, half- 
metallic ferromagnetsii2iiiii5ii& In these materials, one 
spin band is metallic and the other one insulating (100% 
polarized ferromagnet). Since a half metal has a Fermi 
surface only for one of the two spin orientations, it is clear 
that the traditional description for weak ferromagnets is 
inapplicable, and other methods must be employed. 

In Sec. II, we outline the central equations of the qua- 
siclassical theory of superconductivity. Compared with 
the full microscopic theory, the quasiclassical theory of- 
fers considerable simplifications when treating spatially 
inhomogeneous states by reducing the content of (un- 
necessary) information carried by the Green's functions. 
However, this leads to nontrivial boundary conditions 
which have to be formulated at interfaces separating dif- 
ferent metals that connect the solutions of the two sides. 
Such conditions have been derived for nonmagnetic inter- 
faces by Zaitsevfii and for magnetic interfaces by Millis 
ef alm^ After a short description of this work in Sec. Ill, 
we formulate an alternative but equivalent set of bound- 
ary conditions, where the transmission through an inter- 
face is parameterized by a hopping amplitude that con- 
tains the information of various processes contributing 
to particle transfer. This approach enables the formu- 
lation of boundary conditions in a simple and appealing 
form. The equivalence to existing methods is demon- 
strated in Sec. IV. As explained in Sec. V, the advan- 
tages of the t-matrix formulation are especially evident 
in studying interfaces that separate two materials with 
a different structure of the Green's functions and/or a 
different number of trajectories, such as in the case of a 
superconductor/strong ferromagnet interface. Finally, in 
Sec. VI, we apply our theory to study the current through 
a point contact separating a singlet superconductor from 
a strong ferromagnet. 



II. BASIC EQUATIONS OF QUASICLASSICAL 
THEORY 

The quasiclassical theory of superconductivitj*iS*ii is 

formulated in terms of quasiclassical Green's functions 
(or propagators) g(p,R, e,f) that depend on the spatial 
coordinate R and time t and describe quasiparticles with 
energy e (measured from the chemical potential) and the 
momentum direction on the Fermi surface p = p/pf 
moving along classical trajectories with direction given 
by the Fermi velocity Vf{p)}^ The quasiclassical Green's 
functions are 2x2 matrices in Keldysh space (denoted by 
a "check" accent), 
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5^ 

.9^ 



(1) 



with three nonzero elements (retarded g^, advanced g"^, 
and Keldysh g^). In describing superconductivity, these 
elements in turn are 4x4 Nambu-Gor'kov matrices in 



combined particle-hole and spin space (denoted by the 
hat symbol), for example, the retarded Green's function 
has the form 



/ .9ft gfi /n fn \ 

9fi //I //! 

fR fR xR xR 

Jn -'Ti yti 



(2) 



V /it hi 9i\ I 

All these matrix elements are not independent of each 
other. Indeed, the elements in the lower half of the matrix 
are related to the ones in the upper half through the 
conjugation symmetry, e.g. 



~R,A.K I - -D ,\ 



glf'\-p,V.,-e^r. (3) 



The quasiclassical Green's functions satisfy the Eilen- 
berger transport equation 



[ef3-S,5]^+«v/(p)-VR.g = 0. 



(4) 



Generally speaking, the self energy S(p, R, e, t) includes 
molecular fields, the superconducting order parameter 
A = Al, impurity scattering, and external fields. The 
noncommutative product ® combines matrix multiplica- 
tion with a convolution over the internal variables, and 
"^^3 = ■^3! is a Pauli matrix in particle-hole space. The 
quasiclassical Green's functions also satisfy a normaliza- 
tion condition 



(5) 



In addition to Q and ©, self-consistency equations for 
different parts of the self-energy have to be provided; e.g. 
for the (weak-coupling) order parameter the condition 
reads 



A(R,f) = A 



de 



(/^(p,R,e,t))p, 



(6) 



where A is the strength of the pairing interaction, ( )p 
denotes averaging over the Fermi surface, and is 
the particle-hole off-diagonal part of the quasiclassical 
Keldysh Green's function. The cut-off energy Cc is to be 
eliminated in favor of the transition temperature in the 
usual manner. 

When the quasiclassical Green's function has been de- 
termined, physical quantities of interest can be calcu- 
lated; e.g. the expression for the current density adopts 
the form 



^Tr(eiV/v/(p)f35'^(p, R, e, i))p, 



(7) 
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where e is the electron charge and Nf is the density of 
states on the Fermi surface. However, to form a complete 
theory for studying heterostructures, the above equations 
must still be supplemented with the boundary conditions 
connecting the solutions at the separating interfaces. We 
introduce these conditions in the following chapter. 



III. BOUNDARY CONDITIONS 

A. Scattering-matrix approach 

Interfaces represent strong perturbations on an atomic 
length scale and, therefore, fall out of the applicability 
range of quasiclassical theory. However, as was shown 
in the pioneering work of Zaitsev^ interfaces can be 
brought within the quasiclassical theory by means of ef- 
fective boundary conditions that connect trajectories re- 
lated through interface scattering processes. Later these 
conditions were generalized for an arbitrary magnetically 
active interface, i.e. one that scatters quasiparticles dif- 
ferently depending on their spin orientation.^^ The latter 
case is relevant for studying interfaces with spin-polarized 
materials such as ferromagnets. The procedure for the 
derivation of the boundary conditions begins by isolat- 
ing a region of quasiclassical size < 5 around the 
interface located at the origin of the perpendicular co- 
ordinate X {5 much larger than the atomic-size range of 
the strong interface potential but much smaller than the 
superconducting coherence length ^). In the half spaces 
|a;| > (5, the solutions for quasiclassical Green's functions 
can be found by standard methods described in the pre- 
vious chapter. The solutions for the left {I) and right (r) 
sides are then matched via a scattering matrix 
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with cjin = g{p) and (jout = g{p), where p (p) is a unit 
vector along the momentum direction with the perpen- 
dicular component directed towards (away from) the in- 
terface. The boundary condition consists of four cou- 
pled nonlinear equations for the incoming and outgo- 
ing matrix propagators on both sides of the interface. 
Solving this equation system and dealing with the pos- 
sibility of arriving at unphysical solutions is evidently 
not a simple task. Progress towards a more conve- 
nient form of boundary conditions has been made by Es- 
chrig (nonmagnetic interfaces)^ and Fogelstrom (mag- 
netic interfaces). «Si They employed the powerful Riccati 
parameterization method which allows for a considerably 
simpler representation of boundary conditions in terms 
of the Riccati amplitudes>^ However, the conditions in 
Ref. 1^ were only derived for the equilibrium (retarded 
and advanced) propagators. Furthermore, even in equi- 
librium situations they can not be used in the published 
form in the case when the two sides of the interface have 
a different number of trajectories (i.e. when matrices 
Sir and Sri are not invertable). This situation arises in 
the context of half-metallic materials, where trajectories 
exist only for one of the spin orientations. 



the form of which is determined by the detailed micro- 
scopic structure of the interface region and on the qua- 
siclassical level has to be treated as a phenomenological 
parameter of the theory. The crucial simplifying obser- 
vation is that, since the strong (of the order of the Fermi 
energy) interface potential dominates the Hamiltonian 
in the interface region, the scattering matrix corre- 
sponds to that of the normal state, i.e. does not contain 
particle-hole mixing. Also, it has no Keldysh space struc- 
ture. 

The boundary conditions were derived for a smooth 
(on the scale of £,) interface, assuming the conservation 
of momentum p|| parallel to the interface. In the follow- 
ing, all momentum-dependent quantities should be un- 
derstood as having the same p||, unless explicitly stated. 
In terms of quasiclassical Green's functions they adopt 
the form.— 



B. Transfer-matrix approach 

Due to the abovementioned difficulties we proceed in 
an alternative but equivalent routeiS^ This method re- 
quires solving for the auxiliary quasiclassical propaga- 
tors and g^'^ for an impenetrable interface. They are 
to be calculated with the self-energies determined 
with the full propagator, and using the simple perfectly- 
reflecting boundary condition 

gZ = s'gi^{s^)\ (10) 

where i = l,r. They also satisfy the normalization con- 
dition, g*'" ® 5*'° ~ — TT^i. The impenetrable interface 
is characterized by two surface scattering matrices, S*' 
and S"^ . Particle conservation requires them to be uni- 
tary, (5**)^ — (5*)^^. The transmission processes for an 
interface with arbitrary transparency can be taken into 
account with a t-matrix formulation. The transfer ma- 
trices are determined with effective hopping amplitudes 
Tir and Tri by the following equations: 
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= Tlr gout + Tlr gout T/r ® 5i„ ® ^in , (Ha) 
t^n = 5out + ^ri ffo^t ® gin ® > (Hb) 

with Tri — {TirY due to particle conservation. The corre- 
sponding t matrices for outgoing trajectories are related 
to the ones for incoming trajectories through the relation 



matrix in terms of the surface scattering matrices and 
hopping amplitudes: 

\Srl Srr J \ I J \ -r J \ Sr J 

(16) 

where we have defined 



^out ~ ^ '^v 



(12) 



The t matrix describes the modifications of the decou- 
pled quasiclassical propagators due to virtual hopping 
processes to the opposite side. Finally, the boundary 
condition can be expressed in terms of and 5''*^ to read 



gin 

gout 



9r, 



~i,0 
^out 



+ (.9. 
+ {9. 
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tn ' 
out 



i7rl)(g>il^(g>{gl'° -ini), (13a) 
-^^i)®^i„,®(5L°t + ^7^i). (13b) 



In the i-matrix description, the phenomenological param- 
eters containing the microscopic information of the inter- 
face are the surface scattering matrices and the hopping 
amplitudes. The particle-hole structures of the surface 
scattering matrix and the hopping amplitude are con- 
nected through 
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to ensure the conservation of current. In the general case 
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(l+n^firTrl) \l~TT^firfrl), (17a) 

{1 + n'^TriTiry^il - Tr'^TriTir), and (17b) 

{l+Tr'^TirTriy^2nTir. (17c) 



The identity H16(l serves as a precise definition of the aux- 
iliary matrices Si and Sr in terms of the physical parame- 
ters of the full scattering matrix. Using (|15l) , the particle 
{Sp) and hole (Sh) parts of can be seen to be related 
by 



ShiPii) 





Sr 



s;i-P\\) 





Sr 



(18) 



In particular, if the interface scattering matrix is spin- 
inactive, Eqs. (9) reduce to those derived by Zaitsev. In 
the following, we show that the solution of (13) in the ap- 
propriate limit also solve Zaitsev's boundary conditions 
for arbitrary transmission of the interface, and in both 
equilibrium and nonequilibrium situations. On the other 
hand, in the case of diffusive conductors the boundary 
conditions of the t-matrix approach are equivalent to the 
ones derived by Nazarov^ 



^(p,|) = ^*'^(-P||). (15) 

In this formulation, the boundary problem effectively re- 
duces to calculating the auxiliary Green's functions for 
perfectly reflecting interfaces. Numerically this is an ex- 
tremely simple task, e.g.^ employing the procedure of Ric- 
cati parameterization. Afterwards the boundary Green's 
functions for the partially transmitting interface can be 
obtained directly from Eqs. (13), since solving for the 
necessary t matrices (11) only involves a 4x4 matrix in- 
version. When contrasted with solving the group of equa- 
tions (9), the i- matrix approach manifests its usefulness. 



IV. RELATION TO OTHER METHODS 

The underlying perturbative nature of the i-matrix ap- 
proach might arise suspicions concerning its applicabil- 
ity when the interface in question has high transparency. 
The boundary conditions (13) are, however, valid for ar- 
bitrary transmission and, in fact, completely equivalent 
to the corresponding scattering-matrix description (9). 
The connection between the two approaches is estab- 
lished by the following identification of the full scattering 



A. Zaitsev's boundary conditions 

The boundary conditions of Zaitsev read (we suppress 
the symbol (8) and unit matrices for clarity )ii 

gi = 9l = 9a, (19a) 

gaiRigt? + (g;)'] = -^^Dg-gt, (19b) 

where R (D) is the reflection (transmission) coefficient, 
R + D = l, gY = Hg-: - gfut)n, and gf'- = {§1 ± 
gi)/2, with g^-'' = {gl'J^ + g'out)/'^- In the corresponding 
limiting case the surface scattering matrices S^'^ are unit 
matrices, the hopping element can be taken as a real 
number, Tir = t1, and the boundary conditions in the 
t-matrix approach are 

gin = t' + {f + ^^)i'[g'''-^^), (20a) 
glut = + (5'^" - (5'^° + «^), (20b) 

with = gout = g''° and = il^t = i\ The <-matrix 
equations now take the form 
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2 i^l,0~r,0\~1^2~l,0 



From (17) we have 



which, using (18) and the identity 



(21a) 
(21b) 



(22) 



which is identically fulfilled provided that the transmis- 
sion coefficient in the t-matrix description is identified 
as 



l-i? = 



(l+7r2T2)2- 



(30) 



B. Nazarov's boundary conditions 



(23) 



immediately gives H19all . This condition ensures the con- 
servation of current. To show (jl9bp . we first express it 
in terms of the quantities g^^ as follows: 



(1-i?) 



1 + — 9s9s 
ITT I 



, 9a \ ~l 

1 - — I gs9s 



1-1^ 



-2i^{R+l)ga 
where we have used the identity 

i — l,r. Using (17) we find 



= 0, 



~r~l 

9s9s 



whereby Eq. 124|) transforms to 



{(1-i?) 



ga ^ .,0 ...0 

ITT , 



1 + ^ ) 



(24) 



(25) 



(26) 



2iniR+l)ga] ^0. (27) 



This form exhibits directly the unphysical solutions of 
Zaitsev's boundary conditions, determined by vanishing 
of the first square bracket in (|?7|) . The physical solu- 
tions are given by the requirement that the curly bracket 
of H27|l vanishes. On inserting (18) and (19) into this 
expression and using the identity 



The boundary conditions for diffusive conductors, pre- 
sented by Nazarov)^ are formulated in terms of a 
Keldysh-Nambu matrix current, the Keldysh part of 
which defines the electric current through the interface. 
In the i-matrix approach, this matrix is proportional to 
ga of Eq. H22I) and, therefore, to the quantity 



(31) 



determined at the left-hand side of the interface. To sim- 
plify the following expressions, we again choose f;^ — t1 
and real. Furthermore, in the context of diffusive conduc- 
tors both the Green's functions and the hopping elements 
should be regarded as trajectory-averaged quantities, i.e. 
independent of p. Using (21a) and H23|l we obtain 



/ = r'fg\l - r'fg'r' - r'g'g'^a - r^g'^-^ (32) 

where we have dropped the zero from the superscript (all 
Green's functions are auxiliary ones). Writing the matrix 
current in the form 



/ ^ r'fg^i - r'g'gna ~ r^g'^Y^^ ~ r^fg'Y^ 

and exploiting the fact that g'^g^ commutes with 
arrive at 



g^'g'-, we 



/ = -Ag\m-r'g'gl-Hi-r'rgr' 

= -y[g^,g^]{l^T'{g\g'-} + 7rY^)-\ (34) 

Finally, using Eq. H3U|I to identify the transmission coef- 
ficient, and defining = g^'"^ /{in) because of the dif- 
ferent convention for normalizing the Green's functions 
used in Ref. 12J, we arrive at 



(1 - ay^a - (1 - b)-^b = (1 - a)-i - (1 - 6)-\ (28) 
one arrives at the condition 



D[G",G"' 



4. + D{{G\G-'}-2)^ 



(35) 



[(1 - i?)(l + ttV*) - 2{R + 1)^V2] ga = 0, (29) 



which, apart from the prefactor, is the matrix-current 
expression defining the boundary conditions of Nazarov. 



V. INTERFACE PROBLEM WITH 
FERROMAGNETS 

A. Weak and strong ferromagnetism 

As already mentioned, the quasiclassical theory is for- 
mulated in terms of quasiparticles travelling along clas- 
sical trajectories. Smooth interfaces between different 
materials introduce coupling between incoming and out- 
going trajectories with the same momentum parallel to 
the interface. A ferromagnet has a different Fermi surface 
(or, equivalently, set of trajectories) for each of the two 
possible spin orientations. Consequently, two different 
limiting cases that allow a quasiclassical description nat- 
urally emerge (see Fig.nj. In the first case the exchange 
energy splitting of the two Fermi surfaces is small enough 
that the quasiparticle wave packets on the two trajecto- 
ries corresponding to the same parallel momentum but 
different spins overlap and, therefore, the two trajecto- 
ries remain fully coherent in the ferromagnetic region 
(Fig. ^). Technically this means that the full 2x2 spin 
structure of the quasiclassical Green's functions, defined 
by Eq. is to be retained in the ferromagnetic side of 
the interface. This case, relevant for weak ferromagnets, 
has been widely studied in the literature; the standard 
description simply involves a spin-dependent shift in the 
quasiparticle energy, effected by the replacement 

era efs - ha^l (36) 

in the Eilenberger equation Q). Here h is the exchange- 
field parameter and as is a Pauli spin matrix. Other 
Fermi-surface parameters, i.e. Fermi velocities and the 
density of states, are assumed identical for the two spin 
bands in the ferromagnet. 

In this article we restrict ourselves to the opposite lim- 
iting case of strongly ferromagnetic materials, illustrated 
in Fig.^. That is, we assume the exchange splitting and 
the resulting directional deviation of the two spin trajec- 
tories sharing the same parallel momentum to be so large 
that the coherence between them is lost completely. As a 
consequence, the quasiclassical propagators have no ma- 
trix structure in spin space. In particular, conventional 
Andreev reflection processes are forbidden because elec- 
trons and holes in opposite spin bands occupy different 
trajectories which do not interfere with each other. Tra- 
jectories with different spin orientations can only be cou- 
pled incoherently, such as e.g. due to elastic spin-flip 
scattering by magnetic impurities. It should be empha- 
sized that no energy shift of the form (|36|l should be 
introduced in this limit; instead, Fermi velocities and the 
density of states become spin dependent. The reason for 
this is that the integration over the energy of relative 
motion ( integration" ) , employed in the formal process 
of converting the full two-particle Green's function into 
quasiclassical ones, is now performed separately around 
the two different Fermi surfaces. This is in contrast to the 
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FIG. 1: Two quasiclassical pictures of an interface separating 
a spin-unpolarized material (left-hand side of the interface) 
from a ferromagnet (right-hand side): a) weak ferromagnet 
with a small splitting of Fermi surfaces for the two spin orien- 
tations (indicated by solid and dashed curves) and b) strong 
ferromagnet with a large splitting. Incoming trajectory from 
the left-hand side and corresponding outgoing ones with the 
same parallel momentum on both sides are indicated by ar- 
rows. 



case of weak ferromagnets, where the same ^-integration 
range is used for both Fermi surfaces simultaneously. 

A very interesting special case which falls into the lat- 
ter category of ferromagnets with strong spin splitting is 
that of half-metallic materials. In fact, half metals are 
metallic in one of the spin bands only - the other one 
is insulating. Such behaviour has recently been reported 
in CrO-^^'i^ and in certain manganite materials^ and 
has attracted considerable attention because of possible 
applications in the emerging field of spintronics.^^ Since 
in half metals a Fermi surface only exists for one of the 
spin orientations, the standard description for weak ferro- 
magnets is obviously inapplicable. However, half metals 
still allow for a straightforward quasiclassical treatment 
in the separate-band picture: quasiparticle trajectories 
simply exist only for one of the spin orientations. 

B. Spin mixing 

The quasiclassical boundary conditions in the hopping 
description involve surface scattering matrices 5'''' that 
characterize a fully reflecting interface. In the case of 
a magnetically active interface the most general form of 
such matrices (for quasiparticles), satisfying the require- 
ment of unitarity, was pointed out by Tokuyasu et alm^ 
to be 

5 = e-'*/2e-'W2)A-'T^ (37) 

where /t is a unit vector pointing to the direction of the 
surface magnetization and cr is a vector constructed of 
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Pauli spin matrices. The corresponding scattering ma- 
trix for quasiholes follows from Eq. H15|l. Dropping the 
irrelevant overall phase factor $, the surface scattering 
matrix is determined by a single parameter, the spin- 
mixing angle 9. The physics behind spin mixing can be 
visualized as follows: even for a fully reflecting interface, 
incident wave functions penetrate a small distance into 
the forbidden, spin-polarized region. This results in dif- 
ferent matching conditions for waves with opposite spin 
directions and, consequently, different phase shifts for the 
reflected waves. 

The relative phase difference introduced by spin mix- 
ing results in interesting nontrivial phenomena at super- 
conductor/ferromagnet interfaces, even in the absence 
of quantum-mechanical coherence between the two spin 
bands in the ferromagnet. One such example is the re- 
cent prediction of a nonvanishing Josephson current in 
a heterostructure with a mcsoscopic half-metallic piece 
separating two singlet superconductors - driven by spin- 
triplet pairing correlations^ (this effect requires, in ad- 
dition to spin mixing, also the presence of spin-flip cen- 
ters at the interfaces). However, even though spin mixing 
is expected to be an intrinsic feature of any spin-active 
interface, systematic experimental estimations of the typ- 
ical magnitudes of are not yet available. As a guide- 
line for such future experiments, we study in the follow- 
ing chapter the differential conductance of a spin-mixing 
point contact between a singlet superconductor and a 
strong ferromagnet ~ simultaneously offering a view of 
the i-matrix approach at work. 



VI. S/F POINT CONTACT WITH SPIN 
MIXING 



over trajectories with cos 4> > Q. The two spin bands in 
the ferromagnet give two separate contributions to the 
current. From Eqs. (13) follows 



9in 



iJout 



2ni [t, g 



-o-\K 



(39) 



where the t matrix and the auxiliary Green's function g'^ 
(for a perfectly reflecting interface) are to be evaluated 
on the ferromagnetic side where the latter has the simple 



form g^'^ 
with 



-iTT-fa, and g^^° = g^^°F - Fg^'°, 



F 



Fe 
Fh 



tanh I 







2T 



tanh 



V 2T / 



(40) 

where V is the voltage over the contact and T is the 
temperature. We choose the electrical potential to be 
zero on the superconducting side of the interface. Writing 
out the commutator l|39|) . Eq. H38|l reads 



J = E ^ j *(cos0 Tr[i^ - {t^F ~ Ft^)] )%. 

(41) 

Using now Eq. (11), the relation — T^{i^)'^T^, and the 
properties of the trace, we flnd 



j = ^ ireNfv'} I de Im(cos 

Tr[(/^/^)t{)^iV^(F - Fo)] )" , (42) 



We consider a point contact with arbitrary transmis- 
sion separating a conventional singlet superconductor 
and a strong ferromagnetic material. The small (com- 
pared with the coherence length of the superconductor) 
dimensions of the contact and, consequently, the small 
size of the current flowing through it does not apprecia- 
bly affect the state of the coupled half-spaces from that 
corresponding to zero transmission. This offers a simplifi- 
cation by relieving us from the necessity of calculating the 
superconducting order parameter self-consistcntly. Ac- 
cording to Eq. O, the current, calculated at the interface 
on the ferromagnetic side, adopts the form 



E 



9tt)\ 



(38) 



where a =t,i labels the spin band of the ferromagnet, 
each with its own density of states N'^ and the Fermi 
velocity vj. For simplicity, the Fermi surfaces are as- 
sumed cylindrical and the interface specularly reflecting, 
the generalizations are straightforward. The impact an- 
gle 4> determines the angle between the trajectory and the 
current direction. The angular averaging is to be taken 



where we have defined an effective interface potential 
= fg^T^ with g^ = g^^^ the auxiliary Green's 
function on the superconducting side of the interface, 
= (l + iTTfaw^)-!, and Fo = tanh(e/2T). We assume 
that the interface does not flip the spin, i. e. hopping pro- 
cesses from an incoming trajectory in the ferromagnet to 
an outgoing trajectory on the superconducting side are 
without loss of generality determined by two real num- 
bers, fa = Tal, for the two possible spin orientations. In 
this case Eq. (1^ gives 



de 



2ittI Img, 



S 

OLOL 



1 



{Fe - Fo) 



(43) 

where g^^ (gf^) is the 1,1 (2,2) element of the fuU 4x4 
auxiliary Green's function at the interface on the super- 
conducting side. In the presence of spin mixing (de- 
scribed by the spin-mixing angle 9) this has the form 



ecos : 



sin \ 



5tt 



e sin I — cos | 



(44) 
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where fl = — e^, A is the magnitude of the bulk 

order parameter, and g^^ can be obtained by replacing 
9 -9. Inserting into ^3)) we obtain 



J = E ^^/^/ / * ( '^os ^ >+ (^e " Fo), (45) 



with 



|e(sin ■ 



cos f 



) — fi(cos I — in'- 



■T| sm ■ ' 



(46) 

and follows from — > and 9 ^ ^9. For subgap 
energies, |e| < A, jf vanishes because il is real. This sim- 
ply reflects the fact that the contribution from Andreev 
reflection processes vanishes in quasiclassical approxima- 
tion due to the lack of coherence between spin-up and 
spin-down bands on the ferromagnetic side. Introducing 
the normal-state transmission and reflection coefficients 
with Eq. (1201), Eq. can be written for |e| > A as 



2^"V1-(t) 



(47) 



The differential conductance G — dj/dV for \eV\ > A 
can now be obtained by differentiation, and at T = 
adopts the form 
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FIG. 2; The normalized conductance G/Gn as a function of 
el// A for i? = 0.1 (top figure), 7? = 0.5 (middle), and _R = 1 
(bottom). The different curves in each figure correspond, from 
top to bottom, to different values of the spin-mixing angle 
ranging from 6 = Q (dashed curve) to 6 = tt in intervals of 
7r/10. 



G = - ^2e^Nfv'} { coscf) (e = eV) )^ . (48) 

a 

In particular, for a half metal with a conducting spin- 
up band and a reflection coefficient = R independent 
of impact angle 0, the conductance (normalized to the 
normal-state value Gn) reads 

G 



Gn 




(49) 

when \eV\ > A. The contribution due to a finite spin- 
mixing angle 9 has the effect of broadening the conduc- 
tance features near the gap edge. This is demonstrated 
in Fig. [3 which shows the normalized conductance as 
a function of the spin-mixing angle for three different 
reflection coefficients of the contact. In particular, the 



characteristic BCS square-root singularity for a tunnel- 
limit (i? 1) contact is removed. On the other hand, 
for perfectly transmitting interfaces, i? — *■ 0, spin mix- 
ing has no effect. As an additional detail, the maximum 
of Eq. gg, attained at eV/A = {1 + /R)/2i?i/4 when 
6* = 0, is shifted towards higher voltages when 9 > 0, 
vanishing altogether if 9 > Tr/2. 

VII. CONCLUSIONS 

We have presented a quasiclassical theory which is 
suited for detailed studies of heterostructures consisting 
of a wide variety of materials: superconductors (both 
conventional and unconventional), normal metals, and 
both weak and strong ferromagnets. The most crucial 
part of this description is the treatment of boundary con- 
ditions at interfaces separating different materials. These 
conditions are formulated in terms of hopping ampli- 
tudes, containing the information of allowed transmission 
processes, and the corresponding t matrices. Compared 
with the traditional scattering-matrix approach, the t 



9 



matrix approach provides clear advantages for studying 
spin-active interfaces, or interfaces which connect mate- 
rials with different numbers of trajectories or with dif- 
ferent internal structures of their Green's functions. A 
particular example are strong ferromagnets of which the 
half-metallic materials form a special case. In connec- 
tion with such materials, nontrivial physics arises due to 
spin-dependent interfacial scattering processes. The cru- 
cial parameter controlling the details of these effects is 
the degree of spin mixing. At present, there have been 



no attempts to determine experimentally the magnitude 
of this parameter at a spin-active interface. To provide 
a guideline for such studies, and to demonstrate the t- 
matrix approach, we have calculated the differential con- 
ductance for a superconductor/half metal point contact. 
In the tunneling limit of such contacts, the conductance 
depends strongly on spin mixing, and should provide an 
effective means of determining the importance of the new 
physics related to spin-active interfaces. 
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